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ON THE USE OF RESIDUE THEORY FOR TREATING THE SUBSONIC FLOW OF A 

COMPRESSIBLE FLUID 

By Cabl Kaplan ■ 



SUMMARY 

A new mathematical technique, due to Milne-Thomson, 
is used to obtain an improved form of the method of Poggi 
for calculating the effect of compressibility on the subsonic 
flow past an obstacle. By means of this new method, the 
difficult surface integrals of the original Poggi method can 
be replaced by line integrals. These line integrals are 
then solved by the use of residue theory. In this way an 
equation is obtained giving the second-order effect of 
compressibility on the velocity of the fluid. The method is, 
moreover, practicable for obtaining the higher-order effects 
of compressibility on the velocity field. As an illustration 
of the general result, the flow past an elliptic cylinder is 
discussed. 

INTRODUCTION 

There has been in recent years an increasing interest 
in the problems of flow in which the compressibility of 
the fluid is taken into account. The effect of compressi- 
bility on the subsonic flow past various simple shapes 
has been calculated by the method of Janzen and Ray- 
leigh or by a method due to Poggi. Both of these 
methods, however, involve mathematical difficulties, 
which discourage their use for further study of com- 
pressible fluids. For example, the second-order effect of 
compressibility on the velocity field around an elliptic 
cylinder has been calculated approximately by the 
method of Janzen and Rayleigh but, owing to a certain 
limitation in the analysis, the result applies only to 
thick ellipses (reference 1). In order to eliminate this 
weakness in the theory, the calculation was repeated by 
the method of Poggi in reference 2, the result being 
expressed in a closed form. This method, however, 
also involved considerable difficulty owing to the neces- 
sity of evaluating many difficult integrals during the 
course of the analysis. 

In the present paper a method, based on that of Poggi, 
is presented for dealing with problems involving the 
subsonic flow of a compressible fluid. Poggi's method 
consists in regarding a compressible fluid as an incom- 
pressible fluid with a continuous distribution of sinks 
and sources throughout the entire region external to the 
obstacle. In the determination of the second-order 



effect of compressibility on the velocity of the fluid, this 
concept of a compressible fluid leads to a series of double 
integrals extended over the entire region of flow. In 
the present paper, by means of a novel mathematical 
technique utilizing complex notation, the surface inte- 
grals are replaced by line integrals. By means of the 
well-known methods of the calculus of residues it is then 
shown that the second-order effect of compressibility 
on the velocity field around an arbitrary shape can be 
explicitly expressed in terms of residues.. The method, 
moreover, is equally practicable for determining the 
higher-order effects of compressibility. As an example 
of the general result, the flow past an elliptic cylinder is 
discussed. 

It is worthy of mention that a paper which also em- 
ploys complex notation has appeared recently (reference 
3). The purpose of that paper was to complete Hooker's 
treatment of the elliptic cylinder (reference 1). The 
method used is equivalent to the original Janzen- 
Rayleigh process except that the differential equation 
for the velocity potential is expressed in terms of the con- 
jugate complex variables z and z. The treatment in the 
paper is limited, however, to the elliptic cylinder, with 
no attempt made to obtain results applicable to arbi- 
trary shapes. 

THE METHOD OF POGGI 

The equation of continuity for a compressible fluid 
moving irrotationahy in two dimensions and with the 
adiabatic relationship between the pressure and the 
density can be written as follows: 

b<f> bq* .b^ba* 
d'<ft . bx bx by by ,^ 

where 

i> velocity potential 

Co velocity of sound in undisturbed stream 

U velocity of undisturbed stream 

j magnitude of velocity of fluid 

7 ratio of specific heats at constant pressure and 

constant volume 



mm o- *t - i 
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Equation (1) may also be written as 

[' J rKi)!-')> 

where the symbol A denotes the Laplacian operator 

O* o* 
bx> + b? 

and where M is the Mach number Z7/co. It is next 
assumed that <#> can be developed in a series of ascending 
powers of M 2 . Thus 

<t,=4> 0 +4nM 2 +<hM i + ... (3) 

Then if a complex velocity w is defined as 



w=—u+iv= 



it follows that 



bx % by 



w=w a +w l M 2 +w 2 M i + . . . 



(4) 



where' 



and so on. 

bx by 



The magnitude 2 of the velocity of the fluid is given by 

q?=ww=(—u+w)(—u—iv)—u 1 +v 1 
or from equation (4), 

+ (w 0 ^+vbW s +w l Vh)M 4 + ... (5) 

If the expressions for <f> and g 2 from equations (3) and 
(5) are substituted in equation (2) and the terms of 
the same powers of M on both sides of the equation are 
equated, it follows that 



A<fo=0 



(6a) 



. , _lf£>*o b/WoWo\ , b<h d Yl . , . 

^-21"^ bxVW) + ^y by\-TF~)] (6b) 



, lR><fo b/w a Wo\ bfa d/wo^"\l /»„x 

and so on for the higher powers of M. 

It- is to be observed that equation (6a) is Laplace's 
equation for the flow of an incompressible fluid. Thus, 
if the known solution for the incompressible fluid is 
used as the first approximation, the higher approxima- 
tions for the compressible fluid may be obtained suc- 
cessively from equations (6b) , (6c) , etc. Poggi's method 
consists essentially in considering the compressible 
fluid to be an incompressible fluid with a continuous 
distribution of sinks and sources in the entire region 
external to the solid boundary. According to Poggi, 



then, the right-hand sides of equations (6) represent 
successive terms in an infinite series giving this sink- 
source distribution. The first approximation to the 
strength of the sink-source distribution in the plane of 
flow is therefore given by 

If new independent variables z=x+iy and z=x—iy arc 
introduced, then symbolically, 



5a; bz^bz 



~\bz bzj t 



(7) 



and the expression for the strength of the sink-source 
distribution becomes 



kn[bfb%^+WTz^y d y W 



2t 



This expression can be further simplified by means of 
the following considerations: Prom the definition of the 
complex velocity, 



w= 



b<j> ,b$ 



it follows, by use of equations (7), that 
u>=2-^and «?=2-g= 



0) 



In particular, Wo, being the complex velocity of an in- 
compressible fluid, is a function of z only. Expression 
(8) can therefore be written as 



Hz 



)<£r dy 



(10) 



Let the plane z of the obstacle be represented con- 
formally on the plane Z of the corresponding circle. 
Since the strengths of the sink-source distribution 
of corresponding elements of the two planes are equal, 
it follows that the expression for the strength of the 
sink-source distribution of an element of the plane Z 
is given by 



(ii) 



dz dz 



where dx dy has been replaced by — dXdYtsA u\> by 
dZ 

Wq-^t where Wo is the complex velocity of an incom- 
pressible fluid past the circle in the plane Z. 
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If the radius of the circle into which the profile of 
the obstacle is mapped is assumed to be R u the complex 
velocity induced at any point Z P external to the circle 
by a source of unit strength is given by 

1 , 1 , 1 _ 1 , Rf/Zp 



R? 7 ' 
-=-z P 



Z P Z-Z P ^R*- Zp Z 



where there are unit sources at the point Q(Z) and the 

inverse point i?^=-^and a unit sink at the center of the 

circle (fig. 1). It follows, therefore, that the complex 
velocity induced at the point Z P by the sink-source 
distribution given by expression (11) is 




(TP 



Figure 1.— Image of a simple source with regard to a. circle. 

i rrdf s l vz s _ Wt dZdz\ dYdv 




dXdY 



where the integrations are performed over the entire 
region external to the circle. 

As noted before, (Tl T t ) p is the complex velocity in- 
duced at any point Z P external to the circular boundary 
by a sink-source distribution originating in the plane 
z of the obstacle. The actual velocity (vh) P of the 
fluid at the corresponding point z P in the plane z of 
the obstacle is given by 

«*=n^ (is) 

the subscript P having been dropped. It will be observed 
that Wi is expressed as a_ function of the conjugate 
complex variables Z and Z of the plane of the circle. 

THE METHOD OF RESIDUES 

Equation (12) in its present form appears to be 
unmanageable. It can be transformed, however, into 
a form suitable for further treatment by means of the 
following theorem: 




Fiami 2.— Illustration of Stokes* theorem lor the plane. 

Consider a contour in the plane Z of the circle con- 
sisting of a closed curve Ci enclosing the circle of radius 
Ri, a small closed curve d enclosing the point Z P , a 
very large closed curve £ enclosing both Ci and &, 
and the lines AB and CD connecting, respectively, the 
curves C» and Ci, and C 8 and G% (fig. 2). The contour 
is described in such a manner as to leave the area S, 
regarded as enclosed by it, on the left. Now, suppose 
F(X,Y) and 6(X~,Y) to be vector components along the 
X and Y axes, respectively. Then Stokes' theorem for 
the plane states that 

JI(H-oT>* ^^FdX+GdY) 

= f (FdX+ GdY)-^ (FdX+ GdY)-£ (FdX+ GdY) 

(14) 

where the line integrals on the right-hand side are 
taken in the counterclockwise sense. 

If g(Z,Z) is a function of the independent variables 
Z^Z continuous and differentiable in the region S 
enclosed by the contour C, then since 



dZ~bX 



it follows that 



JBy Stokes' theorem, equation (14), it is seen that 
J*JJ|tfX dY=±JjgdY+igdX)=±ijgdZ 



or 



Similarly 

/ l% dX ^-H^+Hf 12 ^^ 2 (i6) 

These equations were first given by L. M. Milne- 
Thomson (reference 4). 
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In the present problem where the area S is the entire 
region external to *the circle of radius B u the curve Ci 
is the circle itself, the curve C± is. a vanishingly small 
circle with the point Z P as center and the curve Ct is an 
infinitely large circle with its center at the origin. By 
means of equations (15) and (16) it is possible to replace 
the surface integrals of equation (12) by line integrals 
around the circles C u C a , and C 3 . For the purpose of 
evaluating the line integrals it will be observed that it is 
possible to make the integrands analytic. The solution of 
any given problem therefore reduces to the evaluation 
of residues, for they form the only contributions to the 
integral of a function that is analytic at all points except 
singularities. In order to transform the surface integrals 
into line integrals, it is first necessary to express the 
integrands as derivatives with regard to either Z or H. 
It is to be noted that the first two integrands of equation 
(12) are already in this form. The last two integrands 
can be brought into the required form in the following 
way: It is observed that 



7 2^5= 



(17) 



The third and the fourth integrands of equation (12) 
become, respectively, 



(18) 



The result of the integrations around the circular 
contours Ci, C», and C» can be explicitly expressed 
by means of Cauchy's theorem on residues. This 
theorem states that, if /(Z) be analytic on a contour C 
and throughout its interior except at a number of poles 
inside the contour, then 



2rxM 



(21) 



where M denotes the sum of the residues olj(Z) at those 
poles which lie within the contour C. Similarly, 



X 



g{Z)dZ= -2riN 



(22) 



By means of equations (15) and (16) it follows that the 
line-integral form of equation (12) is 



where 



(20) 



The integrals of equation (19) are taken successively 
around the circle C t of radius B u the circle C% of vanish- 
ingly small radius B if and the circle C a of infinitely large 
radius i? s . Each of these integrals is taken in the counter- 
clockwise sense, and therefore due regard must be 
given to the sign in accordance with equation (15) j that 
is, the sign of the integrals around Q and Cj is negative 
and that around C s is positive. 



where N denotes the sum of the residues of g (Z) at those 
poles that lie within the contour C. 

Consider, for the moment, onlyjhe first integral of 
equation (19). In general, F(Z,Z) is a function of 
both Z and Z. On the contour of integration C, how- 
ever, Z=iP/Z and therefore F(Z,Z)=F(Z,R I /Z). Thus, 
an analytic function of Z has been created and the 
theorem of residues given by equation (21) can bo ap- 
plied. Similar considerations hold for the second integral 
of equation (19) with Z=jB*/2, but in this case the 
theorem of residues takes the form given by equation 
(22). It is to bo further noted that F&Z) is the con- 
jugate complex of F{ZjZ'). This fact simplifies matters 
to a great extent in that the result of the second inte- 
gration can be written down immediately from the 
result of the first integration. 

Consider the first line integral of equation (19): 
In the integral around the circle Ci_of radius R x with 
center at the origin, the poles of F(Z,Z) associated with 
the terms involving only the variable Z lie within the 
contour, while the poles associated with the terms 
involving only the variable Z lie outside the contour 
because Z=#i s /Z. The poles lying outside the contour 
together with the simple pole at Z=Z P do not contribute 
to the integral around Ci. In addition, since Z=J?j'/Z, 
there will correspond to each zero of 21 a pole of Z at the 
origin. Then, according to equation (21), the result of 
the integral around C t is given by 



r 

Ja 



1 



Z-Zp 



F{Z t Z)dZ=2riS(Z P ) 



(23) 



where S(Z P ) denotes the sum of the residues of 

2^ 2, ^(^'^) at ^ e P 0 * 68 ^''th' 11 tno contour Ci. 

In the integral around the circle Cs of vanishingly 
small radius B t with center at Z=Z P , the only pole 
within the contour is the simple pole at Z=Z P . Since, 
in the limi t as the radius B 2 ->0, Z-*Z P and "2-+Z Pl it 
follows from equation (21) that 



Jo, 



Z-Zi 



■F(Z&)dZ=2TiF{Z P & P ) 



(24) 
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In the calculation of tie integral around the circle d 
of very large radius R s , the first step is to replace E by 
RffZ and thus to render F(Z,Z) analytic. Further- 
more, since the radius R a is ultimately made to approach 
infinity, it is expedient to expand the integrand 

in the neighborhood of infinity. The residue is then 
simply the limit of the coefficient of 1/Z as 
It will be observed that this coefficient is the constant 
term in the expansion of F{Z,R?\Z). This calculation 
can be performed for an arbitrary profile in the z plane 
in the following way: 

. The conformal transformation that converts the 
profile in the z plane into a circle with center at the 
origin of the Z plane and leaves the region at infinity 
unaltered is of the type 



z=Z+a 0 -f-|+gi+ - 



(25) 



where the coefficients do, ck, dt, • . . are, in general, 
complex numbers. 

Now suppose that the undisturbed flow of velocity 
U is inclined at an angle a to the negative direction of 
the real axis and that the circulation r is arbitrary. 
In terms of the complex coordinate Z, with origin at the 
center of the circle of radius R lt the potential function 
of the incompressible flow past the circle is 

/ 0 (Z)=u(ze'-+^§^)+glogZ (26) 
For this flow, the complex velocity is given by 

By means of equations (25) and (27) it can be shown 
very easily that in the limit «> the constant term in 
the expansion of F{Z,R 3 *}Z) is simply XPe*". It follows 
therefore from the residue theorem, equation (21), that 



L 



■ z ± z -F(Z£)dZ=2xiU*e i " 



(28) 



Consider the second line integral of equation (19): 
As noted before, the second integral of equation (19) 
can be obtained immediately from the first. Thus, cor- 
responding to equation (23), the integral around the 
circle Ci is given by 



Jc < Z ~Z~ P 



(29) 



where ^{R*fZ P ) is the conjugate complex of S(Z P ) ex- 
cept that Z P is replaced by R'fZ P and not by_ Z P , and 
where, for convenience, the internal pole at Z=R*[Zp 
has been excluded in fhia calculation. 

The integral around the circle Q yields nothing be- 
cause there are no poles of Z within this contour. 
However, analogous to the residue at the simple pole 
at Z=Zp within the small circle given by equation 
(24), there is a residue at the simple pole at Z=R?jZ P 
within the circle Ci (excluded from equation (29)) 
given by 



j ^^^& Z ) dZ= - 2 ™ F (^ p > Z *) 
Zp 



(30) 



where, in the evaluation of the integral, Z has been re- 
placed by BiVZp and, on account of the substitution 
Z=R l *IZ > Z has been replaced by Z P . 

Finally, for the integral around the circle d of in- 
finitely large radius R s , since P(Z,Z) is the conjugate 
complex of F(Z,Z), it follows, according to equation 
(28), that the residue of 

1 



» Rr 

Z ~Zp 



-J[Z,Z) 



is lPe~ ia and therefore 



J-^ z ^F(Z,Z)dZ= 
c, Z ~Zp 



(31) 



In the summing up of these results, it should be 
remembered that a negative sign is attached to the 
integrals around Ci and Cj and a positive sign to the 
integral around C 3 . It follows that 

where_ the subscript P has been dropped and where 
F(Z,Z) given by first of equations (20) 

f/Ri 7 \ obtained from expression for F(Z,Z) by re- 
\Z^J placing Z by Rf/Z 

S(Z) -given by equation (23) 

obtained from expression for S(Z) by re- 
\ZJ placing Z by B,*/Z 



The last term of equation (32) is equal to — ^ the com- 
plex velocity of the incompressible flow without cir- 
culation past a circular cylinder of radius Ri. 
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In equation (32), the terms in the brackets and the 
last term (representing, respectively, the integrals 
around the circles Ci and C 8 ) are analytic functions of Z 
and are therefore solutions of Laplace's equation. On 
the other hand, the second term in equation (32) 
(representing the integral around the circle C 2 ) is a non- 
analytic function of Z,Z and is a particular solution of 
Poisson's equation. Furthermore, the complementary 
solution and the particular solution are such that, taken 
together, they satisfy the appropriate boundary con- 
ditions of the problem; that is, the velocity of the fluid 
normal to the boundary d is zero and the velocity of 
the fluid at infinity is zero. 

Finally, according to equations (4). and (13), it fol- 
lows that the second approximation to the actual veloc- 
ity of a compressible fluid at any point z external to 
the obstacle is given by 



(33) 



where Wi is given by equation (32) and where W 0 , the 
complex velocity of an incompressible fluid past the 
circular cylinder of radius R lt is given by equation (27). 

In conclusion, it can be stated that, when the set of 
surface integrals given by equation (12) are replaced by 
the solution in the form of equation (32), the labor in- 
volved in the calculation of the second-order effect of 
compressibility has been reduced to a minimum, for it 
is necessary only to evaluate residues. The same proc- 
ess can be used, moreover, to derive equations analo- 
gous to equation (32) for the higher-order effects of 
compressibility on the flow past arbitrary shapes. 

The general result of this paper given by equation (32) 
will now be applied, mainly for the purpose of illustra- 
tion, to the case of an elliptic cylinder. 

APPLICATION TO THE ELLIPTIC CYLINDER 

It is known that the region external to a circle with 
center at the origin of the Z plane is mapped on the 
region external to an ellipse with foci at z=— c and 
z=c by the transformation 



z=Z+^ 



(34) 



Thus, to a circle of radius c/2 with center at the origin 
of the Z plane, there corresponds a line segment extend- 
ing from z=— c to z=e in the z plane; and, to a con- 














0 c.bj 







Figuri 3.— Conforms! mapping of an ellipse Into a circle 

centric circle of radius i?i(>c/2) in the Z plane, there 
corresponds an ellipse in the z plane with foci at z= — c 
and z=c and with semiaxes given by 



c 



(Sec fig. 3.) 



By means of equation (34) the streamlines of the flow 
past the circle can be transformed into those of the 
corresponding flow past the ellipse. The irrotational 
flow of an incompressible fluid past a circle of radius i?: 
is given by the complex potential 

UZ) = u(Ze"+^f^ (35) 

where the undisturbed flow is of velocity U inclined 

at an angle a to the negative direction of the real axis 

and where the circulation is zero. 

Before the various terms of equation (32) for the? case 

•of the elliptic cylinder are computed, the process will bo 

simplified by the introduction of the nondimensional 

It Z z 
quantities ^> and with the letters R x , Z, and 

z being retained to denote the corresponding non- 
dimensional quantities. It follows then from equation 
(34) that 

rfZ 1 
Iz" 



1- 



dZ 



1 — 



z> 



(36) 



and from equation (35) that 
W a =u(e-'°- 



(37) 
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By means of equations (36) and '(37) it follows easily, after the elementary integration indicated in the first of 
equations (20) is performed, that 

^=^-)X^) (1 -^ 

e-^Z+^p+^-RW 1 log §=±j (38) 



-2TP{e ta -R 1 1 e- ta y 



In order to evaluate S (Z P ), given by equation (23), the following expression is needed: 



Z*(Z+ 1) (Z- 1) (Z+Rfl {Z-Rfi 



-2 C7 8 (e^-fi 1 i «-"«)[R 1 *e-^+fi 1 V^+|( e -'--fi I «^)*Z log ^^ ] iZ+1) { Z _ l)t <&) 

In the first term of this expression, there is a double pole at the origin, simple poles at Z=— 1 and Z=l inside the 
circle & of radius B t (>l) , and simple poles at Z= — and Z=Ri* outside the circle d. The poles outside the 
contour Ci do not contribute to the residue S(Z P ). In the second term, there are double poles at Z=— 1 and 
Z=l inside the circle Ci. 
The necessary technique for evaluating the residue at a multiple pole is as follows: 

Suppose the function /(Z) has a pole of order » at Z=a so that (Z— a)"/(Z) is analytic within and on the 
contour enclosing the point a. Then the residue at the multiple pole is 

^J^-^l-. (40) 
By the use of this role, it follows readily from equations (23) and (39) that 

-s^(^^t^Bw{K^^^^-!)-(5^-i-W^i)]) <»> 

where the subscript P has been dropped. 

The expression for SffifZ) follows from the expression for S(Z) if i is replaced by — i and Z by B^/Z. Thus 

The expression for F (R?IZ,Z) is obtained from equation (38) by first forming the conjugate complex 
expression F(Z,Z) and then replacing Z by Rf/Z. Thus, 

(e fa Z*-.B 1 »e- fa ) s , 2R l 2 Z t (Mi*e ia -e- , '')(.e tta Z i +R 1 i e- tta ) 



JrfRi* -\ ™ f V'Z'-RSe-*")* , 2g 1 'ZW g fa - e - fa )(g* fa ; 

U l(Z i -l)(R 1 '-Z*)' i {Rf-Z*)* 



+ 



jj? 1 1 Z 8 (e fa -.B 1 t g- < °)'(g 1 'g fa -g- fa ) T Z-ll (43) 

log^-j 



(fi^-Z 2 ) 2 



* 
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The expression for F(Z,Z) is given directly by equation (38) and can be rewritten in the following 
form: 

_2(a«--ftV-fa)^ s | I p[e-^Z+^^|cB 1 ^-«-*-)' log |=lj (44) 

If the foregoing expressions for S(Z), S{Ri 2 jZ), z), and F{Z,E) are inserted into equation (32), 

it follows that 

I/O* W£ ia \ W,u. lj ft (e ii °^e- tt «)(e { «-R l >e-<°) l ]i , (e t *-R l *e- , «)(Z*e"'+ e -»>) 

--R t Z j_- B? + _ BlF ^_____ _ 

_I(B I ^- e -*)(«'._B 1 V'.)[ log |J±j 

+ 4 ff ' (Z*=S7? lQ g Z=l a — Z^J ^ ~ fi AZ'-1)(Z*-1) : 

+ 2(e'«-i?^-'-)^| 1 p[ <! -^2+^p+|(^'-- c -'-)* log §=f| (45) 

VELOCITY AT THE SURFACE OF THE ELLIPTIC CYLINDER 

The velocity at the circular boundary C x is obtained from equation (45) by putting Z=R i e < * and Z=R,e~". 
Since at the boundary the normal velocity of the fluid is zero, the velocity there must bo purely tangential. 
Then, by means of the relation 

Q u tKV 

where Q u is the tangential velocity at the boundary, it follows from equation (45) that 

fQu\ 1 • ft . x 1» a • (B i o x | p, (g/+l) cos (fl+cfl-ZB.'coB (g-a) . » ■ 4J?. 4 sin* (*+«) * 
\ir) tiret . 2 sm ( fl +«)-2 ffl 6m ( fl + 3 «)+^ i? I *-2B 1 *cos2fl+l an 2a+ a i «-2fl 1 4cos2e+l 

~ 4(^- 2 2f' a cos 2 29+lV (fll4 ~ 1)i ^ ( fl -«)+ 2fl > ! ^ 4 +l) sin 2 cos (0- «)-«,« sin 28 cos (*+«)] logjj£^ 

. 1„ .„ , , , N fii*— 22?! 1 cos 2a+l r/ , „, , . ti* • ,„„ ■ xt i R*+2Ri COS0+1 

+ 4fl.(g. , + D (B ;«_ 2B > cos 2g+l)« [(1 ~ J?I +g '^ 8111 (2»-«)+^ ! 5111 "-■Hi* sul C 2fl +«» lQ g j?>-2Alc03J+l 

+1^.(^-1) ^"^ ^ \^l Y [_fl cos (29-a)+2S 1 I cos a+JtfcOB (2fl+«)] tan" 1 ^§t^ (46) 

As a special case of this equation, suppose a=0; that is, the undisturbed stream is parallel to tin- major 
axis of the ellipse. Then 

(^L^-SW-" *** d - R 'W-» *R7=^R^^ sin 9 

+ J R 1 t sin3flIlog|4±|-2? 1 (5 1 *-l) sin 28 l og ^±^^|±l+2B I [(B I «+l) cos 20-22?,*] tan- ^^"/ 



This expression for the velocity at the circular boundary agrees with the result obtained, in reference 2 (with 
1/2?! replaced by <r). 

Another special case of equation (46) occurs when a—irf2; that is, the undisturbed stream is parallel to 
the minor axis of the ellipse. Then, 

H-^cos 30] log I^-^RjW+I) cos 25-25,1 lo g fr'_2fr 8+ 26 taQ-t ^^] 

, This expression for the velocity at the circular boundary agrees with that obtained in reference 3. 
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According to equation (37) the velocity of an in- 
compressible fluid at the surface of the circular cylinder 
of radius R x is 

f %) =2 sin (fl+'a) 

The velocity at the circular boundary, including the 
second-order effect of compressibility, is therefore given 

by 

(%L.= 2ah < 9 +«>+<%).„„. W 

where is given by equation (46) and where M 

is the Mach number Ufc a . 

Now, according to equations (36), it follows that, for 
the elliptic cylinder, 

dz dZ KzrZr&z* 

On the boundary itself, Z=R\e" and Z=R l e~ a . Then 

It follows, therefore, from equation (33) that, at the 
surface of the elliptic cylinder 



where 



(»). 

m 



(B ] i -2B 1 i cos 28+1) 



is given by equation (47). 



~fu) etTelt 



(48) 



As a final observation it is remarked that the result 
for the example of the elliptic cylinder was obtained by 
the use of the fundamental equation (32) of this paper 
with far less effort than were the results obtained by 
the methods of references 1, 2, and 3. 



Langley Memorial Aeronautical Laboratory, 
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